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The problems of the diffraction of homogeneous plane waves and ground waves by a 
linear shoreline in a planar land-sea surface are discussed. The direction of propagation 
of these incident waves is assumed perpendicular, and that of their magnetic vectors ])arall('l, 
to the shoreline. At the air-land interface, the customary impedance boundary condition 
is imposed while the sea is treated as a perfect conductor; atmospheric and ionospheric 
effects are ignored. Exact integral representations of the solutions are presented. In 
the case of homogeneous plane-wave excitation originating over the sea, the integral repre- 
sentations are employed to obtain expressions for the geometrical optics field and for the 
far-field form of the remaining scattered field, transition regions included. The possibility 
of coastal refraction is discussed. 

1. Introduction 

The first version of this paper appeared in July 1952 as a report [Bazer and Karp, 1952] 
which, however, has been out of print since early 1953. The present version has been solicited 
by the editor, who felt, in view of the continuing interest in the subject matter, that it would 
prove useful to make the essential contents of the report more generally available. At the 
suggestion of the past editor, the authors have extended their original historical survey to 
include mention of subsequent work by other authors on the same problem, and on some closely 
related problems. They have not, however, attempted to improve upon their own results; this 
paper is simply a survey of the results of their report. 

The problems treated arise from the following model for propagation over a land-sea 
surface. The surface of the earth is taken to be the X-Z plane and the region ?/>0, the air 
(see fig. 1). The Z-axis (normal to the j^lane of the page) is taken to be the shorehne and 



PLANE WAVE 
INCIDENT OVER 
SEA 




ay ■ 



AIR (y>0) 




PLANE WAVE Figure 1. Plane-wave excitation originating over 

'^^If""' ^^^'^ the sea, 

LAND 



4^+ah = o 

ay 



the half-planes a;>0, y=^0, and x<iO, y=0, to be the air-land and air-sea interfaces, respectively, 
of the earth-sea surface. It is assumed that this s^^stem is excited by plane waves — ^specifically 
by planar ground waves originating at x= <^ or by homogeneous plane waves originating over 
the land or sea. It is further assumed that these plane waves are polarized with the magnetic 
vector parallel, and with the direction of propagation perpendicular to the shoreline so that 
the total field may be taken to be ''two-dimensional.'^ This field is specified by two vectors, 



1 The research in tliis paper was supported by the Geophysics Research Division of the Air Force Cambridge Research Center, under Contract 
No. AF 19(r22)-42. 
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an electric vector E^==E^(x,?/), parallel to the X-Y plane, and a magnetic vector H = /i(x,?/)Zo,Z3, 
being the unit along the Z-axis. It follows from Maxwell's equations that ^ 

g+P+FA=0, -<.<a.<+<«, 0<y<co, (1.1) 

where 

A-=-[;u«+(4xM<ri/a.)]^/2=|Ar|exp({5), 0<6«1 (1.2) 

c 

and that E^ is related to h by 

E^=-T^r7^77Ti ^oXV,h{x,y). (1.3) 

^aJ[€+(4x(7^/co)J 

In these equations, e and ju are the electric and magnetic inductive capacities; a is the conduc- 
tivity of air; co is the angular frequency; c is the velocity of light; and V^ is the operator 
Xod( )/dx+yob( )/c)y, Xq and yo being unit vectors along the X and Y axes. The total 
field is evidently determined, once h{x,y) is known. 

The sea is treated as a perfect conductor; this assumption leads to the boundary condition 



dh{x,y) 



dy 



=0, x<CO, (1.4) 



Further, if, as is here and hereafter assumed, the modulus of the ratio k/ki of the propagation 
constant k of air to that of the land ki is small and the conduction current is much larger than 
the displacement current in the earth, then the usual conditions on the transition of the tan- 
gential components of the electric and magnetic vectors through a surface of discontinuity 
lead to the following impedance boundary condition at the air-land interface : 



(>Hx,y) 
by 



+ah(x,y) 

yio 



=0, x>0. (1.5) 

yio 
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The quantity a is a complex constant whose argument is t+5, i.e., 

a=|a|exp['i(7r/4+5)]. (1.6) 

In terms of the electromagnetic constants e, ju, a and e^, fii, ci of the air and land respec- 
tively, a is given by 

a=7rrfM~VA~^exp ['i;(7r/4+5)], (1.7) 

where 

rf=c(27rcoM^c70-'/'-X(^/47rVz<^0'^', a;=27rj/. (1.7') 

Here, X is the wavelength in air, d the '^skin depth'' of the land, and v the frequency of 
the excitmg field. If v is less than or equal to 1,000 kc/s and ci is the conductivity of wet 
earth, one can easily verify (cf. Bazer and Karp [1952], Griinberg [1942, 1943]) that the approx- 
imations leading to the impedance boundary condition are valid. Furthermore, it can be 
shown that the ratio 

I (a//:) I < 1/10, (1.8) 

and is 0{v^) for v< 1,000 kc/s. It should be added that the boundary conditions of eqs (1.4) 
and (1.5) are not strictly valid in the neighborhood of the shoreline, so that one should not 
expect the fields derived from h(x,y,) to describe the physical situation in this vicinity. 



2 The time dependence exp (—io^t) and the Giorgi MKS system of units is employed throughout. 
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As indicated above, homogeneous plane and ground-wave excitation will be treated. For 
excitation of the lioinogeneous plane-wave type that originates over the sea, an incident 
magnetic field hl{x,y) of the fonn 

ho(^,y) = ^W [*(-^ cos 6o—y sin do)], 0<^o<^- (1-9) 

is assumed. For excitation of the same character originating over the land, an incident mag- 
netic field of the form 

hi(x,y) = exp [—ik{x cos do+y sin do)], 0<^o<| (1.10) 

is assumed. For excitation of the ground wave type,^ propagating in the negative x-direction 
from x= + co to the shorehne, it is assumed that 

H{x,y) = oxp{-[ix(k'+a')^+ay]}, Im (k'+a')'^>0, (1.11) 

In all cases the corresponding electrical excitation is to be determined from eq (1.3). 

To complete the formulation of the problem, it is necessary to specify appropriate condi- 
tions at the shoreline and at infinity. At the shoreline one assumes, whatever the excitation, 
tliat the total magnetic field h(x,y) satisfies 

Lim fVUcos^, psiu^) ^^(^'•»^^- P^"^^) prfe^o. (1.12) 

p->0 Jo <JP 

Here, h* denotes the complex conjugate of h, and a^=p cos 6, y=p sin d,0<6<7r. This condition 
is equivalent to the statement that electromagnetic power is neither radiated nor absorbed at 
the shoreline. 

Since it has been supposed that the air is slightly conducting (see eq (1.2)), it is to be 
expected, on physical grounds, that no scattered power will reach infinity. Writing 

h(x,y)=ho{x,y)+v(x,y), (1.13) 

where hQ{x,y) is any one of the incident waves described above, one finds easil}^ that no 
scattered power will j'each infinity if the following requirement is met: 

T • r^'' */ n ' M c)i<p COS ^, p sin ^) ,^ „ /I i^N 

Lim ?;*(p cos 6, p sm 6) — — ^-^-^- pdd=0. (1.14) 

p-^oo Jo ' ^P 

The problem is therefore to find solutions h(x,y) of the form given in eq (1.13), which cor- 
respond to tlie various excitations introduced above and which satisfy: (1) the time-reduced 
wave equation (1.1), (2) the boundary conditions of eqs (1.4) and (1.5), and (3) the conditions at 
the shoreline (eq 1.12) and infinity (eq 1.14). It has been shown by Bazer and Karp [1952] 
that such solutions are uniquely determined. Here, it may be mentioned that neither existence 
of the solutions, nor questions of uniqueness depend critically upon the fact that arg (a) = 
(7r/4) + 5. For the unique existence of the solutions, it is suflScient to require merely that arg (a) 
be restricted to the range 5<Carg {a)<C.Tr-\-8, Similarly, the small positive imaginary part of 
\k\, k, namely |A:| sin 8 (cf. eq (1.2)), may be dispensed with and (1.14) replaced by the Sommei- 
feld condition. It is, however, retained here because its presence shortens the anal.ysis. 

For plane wave excitation originating over the sea, the formulation sketched above is a 
variant of that given by G. A. Griinberg [1942, 1943]. In these works he considers only the 

vertical component of the electric field, Ey, which from eq (1.3) is proportional to our ^j 

3 The concern licrc is only with the mathematical possibility of such excitation, not with its physical realization. 
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Ey clearly satisfies the wave equation and the same boundar^^ conditions as h. Griinberg 
reduces the problem of finding Ey to that of solving an integral equation of Wiener-Hopf type. 
While recognizing the possibility of exact solution by the Wiener-Hopf technique, he prefers 
an approximate method which leads directly to an expression for the far field of Ey on the land 
and to the conclusion that there is no refraction at large distances from the shoreline. 

Problems closely related to those formulated above have been treated by several investi- 
gators, all employing the Wiener-Hopf technique. Among these are works by S. Edelberg 
[1952], T. B. A. Senior [1952], A. E. Heins and H. Feshbach [1954], and A. F. Kay [1957]. 
Edelberg^s work relates in part to the diffraction of a normally-incident plane electromagnetic 
wave by an imperfectly conducting half -plane. Senior treats the same two-dimensional 
problem but allows arbitrary angles of incidence. Kay's investigation is a generalization of 
the present one in the sense that he employs the impedance boundary condition bh/c)y-\-aQh=0, 
aQT^O when x<iO and c)h/dy-\-aih=0, ai9^0, aiv^a^ when a;>0. Kay, however, assumes that 
ao and ai are real, whereas here the analog of ai, namely a, is complex and ^o vanishes. Heins 
and Feshback [1954] discuss the coupling of two half-planes. The common feature of all of 
these works is the '^factorization^' of a function (j{v) the form a-(j/)=l + K/yF— j'^ (see section 2) 
where i' is a complex variable and k is a suitable constant. 

Concerning treatments of the land-sea problem not using the Wiener-Hopf technique, 
mention may be made of the works of P. C. Clemmow [1953], H. Bremmer [1954], K. Furutsu 
[1955], J. R. Wait [1956, 1957a, b, c, 1958], E. L. Feinberg [1959], and T. B. A. Senior [1956]. 
The reader is referred to these papers for additional bibliographical material and for the latest 
developments in the subject. 

In the next section, solutions of the problems formulated above are given. These solutions, 
which are expressed in the form of integral representations, were originally derived by Bazer 
and Karp [1952] by a variant of the Wiener-Hopf procedure (see, in this connection, the work of 
S. N. Karp [1950a, b]). Here, for the sake of brevity, only the end results of the procedure, 
the integral representations, are presented. It is then verified by standard function-theoretic 
techniques that the integral representations actually do furnish solutions of the problems con- 
sidered. An excellent survey of methods based on the Wiener-Hopf technique for the solution 
partial differential equations has been given by B. Noble [1958]. 

In section 3, the last section, formulas for the far field resulting from plane-wave excitation 
over the sea are summarized. In the special case of horizontal incidence, the far field on the 
air-land interface is shown to agree with that obtained by Griinberg in the works cited above. 
It is noted, in addition, that the behavior of far fields on the land for intermediate distances 
indicates the possibility of coastal refraction. All comments on this subject are, however, 
necessarily of a tentative nature, since, among other things, the vector character of the fields 
is not taken into account in the present two-dimensional treatment. For a thorough discus- 
sion of coastal refraction and related phenomena, based on methods which do not make explicit 
use of the Wiener-Hopf technique, the reader should consult the papers of P. C Clemmow, 
E. L. Feinberg, and T. B. A. Senior, that were mentioned above. 

2. Integral Representations of the Solutions — the Factorization 

2.1. The Integral Representations 

Hereafter, the symbols P^ and Pi will be used for abbreviated reference to the problems 
associated with homogeneous plane-wave excitation originating over the sea and land, respec- 
tively (cf. eqs (1.9) and (1.10)), and Pg will be similarly employed for the problem associated 
with ground-wave incidence. Let m, r, t, and h be defined by the equations 

7n=k cos 6*0, 0<^o<7r/2 (2.1) 

^ Jk.me^-a ^ (2.2) 

%k sin 6i)-\-a 
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<=2-H-(l-r)(F-m2)i/2=J^iyiM_; (2.3) 

ik sm ^o+a 

6=(F+a2)>/2^ Im(F+a2)"2>0. (2.4) 

Let, in addition, (x*{v) and (t~(j') denote the "factors" of the function 

<j{v) = l+[ali{F-vy\ (2.5) 



in the sense that 



= <v). (2.6) 



Explicit representations of these factors will be given later. For the following, it need only be 
required that o-'^iy) and (t~{v) enjoy the following properties: (1) (t^{v) is regular, zeroless, and 
uniformly bounded in the (upper) half-plane defined by Im {v)^—\k\ sin 5; (2) (j~{v) is regular, 
zeroless and uniformly bounded in the (lower) half-plane defined by Im (^)<C \k\ sin 5. Then, in 
terms of these quantities, the integral representations of the solutions may be expressed as 
follows : 

P: :h^=2oxp (imx) cos [^(F-m^)i] .-^y-. C ^^^f^^^r^^^^ (2.7)- 

7^^o-+(m)J_„ {v—m)[i(k^—vy] v/ ^ v / 

P+ :/i^=2exp (iinx) cos [7/(F-m^)^] .-^ C ^fP^'^^f'^^^f'^'^l^ll <T+(v)(h; (2.7) + 

Pf :/i^=exp (— 'imx){exp [— 'i?/(/:^— m^)^]+r exp [^^/(F— m^)^]} 

< f- exp{i[.:B+y(F-.^)J]} 
~:ricr-(-m)J_„ (. + m)[i(F-.^)i] <^ W^'", (2.S) 

P(+ : A'=exp (—imx) {exp [— ii/(F— TO^)i]+r exp [i2/(F— m^)*]} 

< f" exp{i[.x+?/(F-.^)i]} + , , , ... ^- + 
~rt^-(-m)J_ (.+ m)[a+i(/t^-.^)i] '' W''". (2-«) 

P, :A^=exp [-(^6a;+ay)]+-^^J_^ (^^5)(p_y2)i '^ W^/"; (2.9)" 

P; : A^=exp [-i^bx+ay)\—^j_^ (.+ 6)[a+i(F-.0^] '-"('-)^^''- ^2.9) + 



In these equations, the P"^ representations are obtained from the corresponding P~ repre- 
sentations simply by replacing a~{v) by the equivalent expression (T'^(v)[l-\-a/i(k^—v^)^]~^ (see 
eqs (2.5) and (2.6)). In each equation, the terms in the right member consist of an incident 
wave and when reflection is possible (as is the case in the Ps and Pi representations) a suitable 
reflected wave. Thus, in eq (2.7), the term 2 exp (imx) cos [y(k^—m^f] is a superposition of the 
incident wave ^0= exp [ik(x cos Oo—y sin do)] (cf. eq (1.9)) and the wave exp [ik{x cos ^o+^sin^oL 
reflected off the perfectly conducting sea. In eq (2.8) the first term is the incident wave 
originating over the land and the second term is the wave reflected off the imperfectly conducting 
air-land interface, the coefficient r being the reflection coefficient at this boundary. The 
exponential terms in the integrands of eqs (2.7)-(2.9) are factors of product solution, obtained 
from the wave equation by means of the standard separation-of- variables procedure. The 
remaining factors in the integrand are chosen so as to meet the boundary conditions at tlie air- 
sea and air-land interface (see eqs (1.4) and (1.5)). 

To insure convergence of the integral representations and at the same time to satisfy the 
^'damped radiation condition" at infinity (see eq (1.14)), it is necessary to fix the branch in such 
a manner that exp [iy{k^—v^f] is exponentially damped whenever v is on the contour, 
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— oo<^Re v<^co ^ and \v\ exceeds \k\. For this purpose, one fixes the value of (k^—v^y at 
V =iQ hy requiring that 

{k'-vy\.^o=+k. (2.10) 

Choosing the branch cuts along the rays which go from +fc to «> exp (i8) and —kto — «> exp (i5) , 
one may then write 

v—k=ri exp [i(yi+d)], 0<7i<27r; 



where 



v+k=r2 exp [i{y2+d)], — 7r<72<7r, 

ri=\v-k\, 7i=arg (v— t)— arg 5, 
72=arg (z^+^) — arg 8. 



r2 



\p+k\, 

From eq (2.11), it then follows that 

{k'-vy^=^{nr2)^ exp [i(Ti+T2+37r/2+6)], 



0<Ti<27r, 
-7r<^72<7r. 



(2.11) 



(2.12) 



(2.13) 



It can now be verified that arg (k^—v'^)^^Tr/2 for v real and |i^|>^. The behavior of arg (k^—v^)^ 
on the Re (j/)-axis and immediately above and below the cuts is depicted in figvire 2. 

It has thus far been assumed in the Ps and Pi representations that dQ7^T/2, Physically, 
this angle 6q='kI2 corresponds to plane-wave excitation incident normally on the land-sea 
surface. This restriction may be removed as follows. Starting with the P^-representation of 
eq (2.7)"^, one allows ^o to increase to 7r/2. The pole in the integrand at v=m=k cos ^o then 
approaches to contour, on the Re v axis from above. In the limit as is well-known, half the 
residue of the pole at i^=lim m = is split off and the resulting integral is taken as a principal- 

eo-^ir/2 

value integral with respect to the pole at v=0. Letting h^(x,y) be the limiting solution, one 
obtains the following integral representation : 



P+ : h''=2 cos ky- 



a exp (iky) 



{a+ik) Trio-" 



^cr+(0)J_, 



exp {i[vx+y{k^—v^)^] 



(T+{v)dv, (2.14)^ 



the principal value being taken with respect to the pole at zero. P^ refers to the problem 
associated with normal incidence and the function a'^iy). Writing (t'^{v) in terms of (j~{v) by 
means of eqs (2.5) and (2.6), one then obtains the alternative P~ -representation; namely. 



Pn'.h''=2 cos ky- 



a+ik 



exp {iky)- 



^+(0) j-co 



exp i[vxAry{¥ 



v[i{¥-vY\ 



^a-{v)dv. 



(2.14)- 



It can be shown, using a similar procedure, that the same results are obtained when the P^- 
representation is taken as the point of departure. 

It is now a simple matter to verify that the integrals appearing in eqs (2.7), (2.8), (2.9), 
and (2.14) are properly convergent and admit all operations under the integral signs, neces- 
sary for proving that the boundary conditions and the wave equation are satisfied. First, 
taking into account the boundedness of |c7+(i')| and |o-~(v)| on the contour (see properties (1) 
and (2) in the neighborhood of eq (2.6)), the behavior of {¥~v^)^ on the contour, one sees that 



^/iii~i^^Tr 



A 



y^. 




.^^ 




Figure 2. The cut v-j)lane. 
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the intoonuids are at worst 0(v~^) as |j;[— >oo^ independently of x and y. It follows that all the 
integrals eonv^erge uniformly for all x and y in the ranges —^ <Cx<C^ , <^<C°° *^^^d that taking 
(he limit y—>0 under the integral signs is permissible operation. Next, when y is positive, the 
function exp [iy{k^—v^)^ decays exponentially as |?/|-^oo on the contour. This implies that 
one may differentiate any number of times imder the integrals signs and still retain the ab- 
solute convergence of these integrals for all x and ?/, — oo <^x<^oo ^ 0<C2/<C°° ; the reduced-wave 
equation is therefore seen to be satisfied. Finally, if in P'^(P~) representations, all contours 
are bent upwards (downward) slightly so that, as li^l-^oo on the contour, the relation 0<CTni v<C, 
\k\ sin 5(0>Im v^\k\ sin 6) is satisfied, then for xy>0{x<CO) the function exp (ipx) is exponen- 
tially damped. It follows, in all cases, that diflFerentiation under the integral signs is a per- 
missible operation for all ^^0 provided |:r|>0. Employing the P^ represeiitations when 
xy>0 and the P~ representations when x<CO, it is now a routine matter to verify i\\o l)0undary 
conditions. 

2.2. The Factorization 

The crux of any boundary-value problem which is solved b}^ the Wiener-Hopf technique 
is, as a rule, the factorization of a suitable function — aiv) in the present case. In fact, provided 
that the factorization is feasible, one can learn with a little practice how to write down integral 
representations of the solution (e.g., those of subsection 2.1) at once without resorting to the 
usual analytic procedures encountered in the early applications of the Wiener-Hopf method. 
Such an approach is exemplified in the papers by J. Bazer and S. N. Karp [195G] and S. N. 
Karp [1957]. Now as is well-known, formulas for the factors of a rather large class of functions 
have been given by N. Wiener and E. Hopf [1931]. Tlie existence of factors with the desired 
properties is thus assured. However, these fornmlas express the factors in terms of rather 
complicated integrals involving contours of infinite extent and it is often found, for this reason, 
that the behavior of the factors in regions of interest is obscured. Thus, althougli possil)ility 
of factorization is in principle guaranteed, the problem of transforming the original WicMu^r- 
Hopf formulas into a more tractable form remains. This transformation was descrilx'd in 
detail by Bazer and Karp [1952] for the function 

cT{v) = l + [afi{k'-vY^]. (2.5)' 

Here, only the end results of this transformation will be presented and some direct consequences 
of these results will be summarized. 
First, it is to be noted that 

<^~W=^ (2.15) 

can be shown to hold for all v so that it is oniough to give the form of a'^{v). For this purpose, 
let the functions/(c) and L'iv) be defined by 

... 1 Ft . ,,,,-] 1 fTT, ., riv±(^-f}n\ arc cos ju/k) 



and 



(2.16) 

x.(.)^=g_,r .f(-)-f) +. /(^)-/(-^) i (2.17) 

v+k L ^ — ^ v-to J 

where 

6-(F+a2)^, Im 6>0. (2.18) 

In"eq (2.16), the principal values of the inverse trigonometric function and of the logarithm are 
intended and the branch of (F— r)^ is specified as in eq (2.10). In terms of these functions, 
a'^iv) may be expressed as 

(T+(.)-cT+(0) exp J^^. j^U{t)d^y, (2.19) 

wliere 

a+iO) ==[i+a/ik]i. (2.20) 
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Since /(j') is regular at all points of the z^-plane except at the lower branch cut of {¥—v^)^ [cf. 
figure 2], it follows, in virtue of eq (2.17) that L' {v) is regular in this region. The same is 
therefore true of (7'^{v) which may now be shown to have all the properties required of it in 
subsection 2.1. 

The formulas for the far fields of problem Ps, Pi, and P^ involve the values of (7+(fc cos B) 
as 6 varies in the interval <^<7r. It is therefore desirable to obtain a reliable approximation 
for (7+ for d in this range. That it is possible to do so is a consequence of the fact that L'(f) = 
U[^,{alk)] (see eqs (2.16)-(2.18)) is a regular function of ajk when —l<i^lk<l and the fact 
that \alk\ <0.1 (see eq (1.8)). Expanding the exponential term in eq (2.19) as a power series in 
ajk, one finds that 



(T 



■Ik) 



where R'^{9,a/k) is the remainder of the series within the braces. In the same way, it is found 
that 

a-(k cos ^-«/^)=-?(^{i+[^(i-ii^)] w^)-^ [Ki~ii^y 

~i coT(x/2) ]^'^^^^'+^"^^'"^/^^}' X=x-e,0<X<x. (2.22) 

For 6 and x in the indicated range greater accuracy could be achieved by including more terms 
of the power series. However, in deriving the air-land boundary conditions relative errors in 
the fields of the order of magnitude of 0.01 are neglected. Greater accuracy would therefore 
be physically meaningless. 

Bazer and Karp [1952] have shown that 

l.+iMMl^lMm^om, o<9<uo°, (2.23) 



I (7+ (A: cos e,a/k)\ 
\R-(x,a/k)\ 



\(T+(P,a/k)\\a-i—k cosX,a/k)\ 



<0.0], 0<x<140°. (2.24) 



From the relation a'^(k cos e) = {a-\-ik sin 6)a~{k cos d)/i sin d (cf. eq (2.5)), it follows easily 
that the relative errors in 0-+ and a~ when calculated by means of eq (2.28) and (2.24) are less 
than 0.01 throughout the angular range 0<C^<C7r. 

Formulas for (t^(v) of the type given in (2.19) have also been derived by T. B. A. Senior 
[1952] and A. E. Heins and H. Feshbach [1954]. A derivation is also given in B. Noble's book 
[1958], pp. 91-92. It does not appear to be generally known, however, that V. Fock had 
already given the following compact formula for a'^(v) as early as 1944 [Fock, 1944, p. 45]: 



_L/7 ^N /cos ^+cos a ^^.^ r 1 r^+" u , "1 /n ^rN 

(T+{k cos 6)=^ —r-7-^ — — - exp^^ — -, du >- (2.25) 

^ V 1+cos^ \^ ^TjB-a smu J ^ 

Here, the angles and a are defined by 

v=k cos 6, (2.26) 

sin a=ia!k. (2.27) 

Fock merely states the result; but it is not difficult to show that this expression follows directly 
from eqs (2.19) and (2.17) on substituting cos 6 for v/k and sin a for ia/k (cf. eqs (2.26) and 
(2.27)). 
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It should be stressed finally that, in the problems under discussion, the constant a is 
assumed to have nonvanishing imaginary part (see eq (1.6)). If a is real then the absolute 
values of a'^ik cos 6) and a~{k cos d), 0<^<7r, reduce to simple expressions. For, in this case, 
the angle a in eq (2.27) is purely imaginary and 6— a, the lower limit of integration in eq (2.25), 
is the complex conjugate of 6-\-a in the upper. It follows easily from this fact and eqs (2.15), 
(2.25)-(2.27) that 



\(7-{k cos6>)|- 



I (7+ (A: cos (jr—d)] 



cos ^—1 



COS ^— COS a 



Iv—^k'+a' 



(2.28) 



whenever v/k is real and \v/k\<l. These formulas were first derived by A. F. Kay [1957, eq 
66]. Kay incorrectly implies that this shnplification is valid for nonreal values of a, specifically, 
that eq (2.28) is also applicable to the problems discussed hei'e. 

3. The Far Fields-Coastal Refraction 

3.1. The Far Fields 

For the sake of brevity, only problem P^ will be discussed. Tlie corresponding results 
for the remaining problems are given by Bazer and Karp in their report [1952]. 
Write 



x=pcosd, y=p sin 6, 0<d<T, 
and introduce the new variable of integration 

i3=arc cos (v/k), 0<i3<7r 
into eq (2.7)^. The result is 
Pt :A'=exp [ikp cos (<9— 6'o)]+exp [ikp cos (d+do)] 

exp [ikp cos {fi—d)]a~^{k cos (3) sin fidp 



(3.1) 
(3.2) 



__±_c 

TTia-^ irn) J c 



(cos /3— cos ^o) ia-^-ik sin jS) 



(3.3) 



wliere the contour (J is the image in the /:^-plane of the contour Im ^^=0 in the ^-plane. It is 
easy to verify that C is essentiall}" the contour shown in figure 3. C is symmetric with respect 
to the point ,3=7r/2, and as Im fi-^i^ on C, Re p-^b=^vgk. The lower portion of the upper 
cut in the i/-plane goes over into the upper half of the line Re 1^=0 (darkened in the figure), 
whereas the upper portion of the upper cut goes into the lower half of the line Re fi=^. The 
segment traced from to +^ in the j'-plane goes over into the arc traced from 7r/2 to in the 




FiGUBE 3. The complex ^-plane. 



fl=(7r-S)-ico ^=(e + |;)_|co 
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/3-plane. Symmetry with respect to /3=7r/2 determines the rest of the figure. The deviation 
of C from the darkened rectihnear path, which is due to the presence of a small imaginary 
part in k, has been greatly exaggerated for the sake of clarity. The poles in the integrand 
of the P^ representation are at 13=00 <t/2 and at /3&=arc sin (ia/k)^ (ia/k) = \a/k\ exp (i3w/4:). 
These poles are shown enclosed in full circles. One might expect a pole at P'^ir—iialk). 
However, on substituting v=k cos ^ into eqs (2.5) and (2.6) it is found that 
i sin jS (7+ (A: cos /?)/ (a-\-ik sin fi) = a~{k cos /3). Thus the integrand in the P^ representation is 
regular near /3=7r (since a~{k cos /3) has this property). Note, since arg (ia/k)^37r/4, /?& is out- 
side the strip < Re jS < tt. Nevertheless, this pole makes its presence felt owing to its prox- 
imity to the saddle point at ^=0.^ 

Now that the positions relative to C of the poles at l3=do and (3=?irc sin (ia/k) have been 
established, it will be supposed, for the sake of simplicity, that 5= arg k is equal to zero. The 
path C then becomes the heavy rectilinear contour of figure 3 except at the pole |S==^o which 
it avoids by means, say, of a small semicircular arc. The shaded portions of the figure are 
those regions where, for fixed 6, Im [i cos {^—d)]<CO; one may therefore deform the contour 
freely in this region. 

The steepest descent contour S{6) (see fig. 3) is defined by the equation 

S{6) : cos (13-6) = 1 + is', (3.4) 

where s is real. Solving for 5 in terms of /3 one finds on selecting the proper branch that 

s=+2'^ exp (iir/i) sin f(0-^)/2]. (3.5) 

By deforming (7 into S{d) one obtains from eq (3.3), in addition to the appropriate com- 
bination of incident and reflected waves, an integral term of the form 



J s 



exp [ikp cos (l3-e)]A{^)dl3 (3.6) 

Sid) 

which, in virtue of eq (3.4), becomes 

J{e)=2i exp [iikp-w/^)] f exp {-kps')B(s)ds (3.7) 



where 



_^_[^(£)l 



Let k be assumed fixed and let tp>>l when p>Po; here po is a sufficiently large positive 
number. Then, in the evaluation of the integral term of eq (3.7) only the values of B(s) in 
the neighborhood of s=0 are important owing to the presence of the exponential factor 
exp {—kps^). To obtam a reliable and tractable approximate expression for J (6), one there- 
fore is led to replace B(s) by an approximation chosen so as (1) to represent B(s) well near the 
origin and (2) to make the resulting term wise integration simply evaluable in terms of known 
functions. Suppose first that B(s) is regular in a sufficiently large neighborhood of the origin. 
In this case a power series expansion for B(s) meets the above requirements, since it leads to 
integrals of the form 

exp {—kps^)s^''ds 



r. 



which can be evaluated in terms of inverse powers of kp with gamma-function coefficients. 
The leading term of this expansion gives the same result as the usual saddle point approxima- 
tion [s=0 corresponds to the saddle point at 13=6^ since 5=(l+'i) sin |(/?— ^)]. On the other 
hand, suppose B{s) possesses a pole at s=Sq, where \so\ may be made arbitrarily close (or even 



^ It should be recalled here that \alk\ <1/10 in the range of frequencies studied (see eq (1.8)). 
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equal) to zero. Here, for fixed kp>kpQ^y>l a series expansion may possess too small a radius 
of convergence to extend readily into the range where the influence of the exponential decay 
of exp (—kps^) is felt. In this case, following the methods of H. Ott [1943] and B. L. Van der 
Waerden [1951], one isolates the pole in a suitable fashion. The method of Ott differs from 
that of Van der Waerden in the manner of isolating the pole. In their report [1952], the authors 
employ the procedure of Van der Waerden. The following is a summary of the end results 
of this procedure for problem P^. 

3.2. Summary of Results 
a. Characterization of the Regions of Figures 4(a) and 4(b) 

The poles at l3=do and P=fib are the poles which play a role in the far field formulas. A 
careful examination of the relative positions of S(d) and ^^ reveals that S{d) sweeps over 
jSft as d approaches zero.Z The value of 6 for which S{d) actually passes through /?;,, is called On. 
Expressed in terms of the variable s, the poles at 13= 6[ and /?=/?/, are located respectively at 
(see eq (3.5)). 



and 



^0=2* cxp (i7r/4) sin [{do-d)/2] 
So=2i exp (7r/4) sin [(P,-e)/2], 



(3.9) 
(3.10) 



Let tp>frpo>>l. Then the value exponential term is appreciable for those values of 6 such 
that /:p,s'o<C<Cl and the Van der Waerden procedure applies. Let it be assumed for the present 
that 007^0. The region included within the parabola (see figure 4(a)) 



B: 



\{2kpy^sm[{do-d)/2]\=K,<<l 



is called region B, Similarl}^, the equation 

D: \{2kpy^ sin [(/?,- ^) /2] |= if 2«1 



(3.11) 



(3.12) 



defines a parabola-like region D in the vicinity of B=Bi,. In figure (4a) the region D is indicated 
roughly by the shading above the positive a:;-axis. Since ^d has a small imaginary part, this 
region, unlike region B, cannot extend to infinity. The regions A and C of figure (4a) are 
defined bv the relations 



A: 



^>0>Oo+Ado 



(3.13) 
(3.14) 



where A^o <^iid A^^, are small positive angles. Since (3 is bounded away from the poles when 6 
lies in A or C it is possible to choose po so large that \(2kp) sin [(^— ^o)/2]|>>l and \(2kp)' 
sin [((Sft— 0)/2]|>>l whenever p>po. In these regions the usual saddle point approximation 
applies. 



DIRECTION OF 
REFLECTED RAYS 




Figure 4(a) Regions A, B, C, D — characterized by 
differing asymptotic field behaviors. 

The excitation is assumed to originate over the sea and to make an 
angle of ^o, O<0o<7r/2 with the positive j-axis. 



a This fact is proved as follows. Let fib=Xb-\-Wh xi, t/b real and w^rite ^=i:o+?>o for that point on S{0) which has the same ordinate as /Si, so 
that y^=yh. It is then easy to show that x^<.Xb. 
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DIRECTION OF 
INCIDENT WAVE 






„,|..„i.„MH, 






B':j(2k^)%in(-|)|<K3 



Figure 4(b) Regions A' and B' — characterized by 
different asymptotic field behaviors. 

The excitation progresses parallel to the sea towards the shoreline. 



Figure (4b) is obtained from figure (4a) by setting ^o=0- Regions A^ and B^ are defined 
by eqs (3.13) and (3.14) with i9o=0. 

b. Notation 

The following notation will be employed in the next two sections. 



2i exp (iT/4:){a—ik sin fi)a+{Fcos jg) sin cos [(fi+do—2d)/4:] 
2 cos [(/3-^o)/4] sin [{P+do)/2] cos [{0-do)/2] 

h, = h,{e,6o)=2^ exp (i7r/4) sin [%-e)l2l 

l\=l\{efi,) =2^ exp (i7r/4) I sin [{6,-6) I2]\) 

rr... , 2^ exp (i7r/4) cos^ (g/2)(7+(/: cos /3) sin /3 cos [(g-2^+fe)/4] 
^2lP,t/,(^o; ^^ ^^g [(^-^,)/4)] cos[(^+^,)/2] cos![(^-^)/2] 

b2=b2{6M=2' exp (^7^/4) sin [{p,-6)/2l 



sm I3^=ia/k, cos 18^,=^ — j — -) 



11=11(6 A)=T2'^ exp (W4) sin [(/3.-^)/2],| '^^^'^ ^>^^'*^ 



Im cos i5&>0, 

— when ^>^ 
+ when 6<C6b]^^ 



1(: 






>o 
<o 



The constants c and r are defined by the equations 

c=a/{a'+k' sin 6) 



(3.15) 

(3.16) 
(3.17) 

(3.18) 

(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 



r= (ik sin 6o—a)/{ik sin ^0+^). 
The function erfc (z'), with z complex, is defined by 



erfc (^^) = l-erf {zh)=l^T,-h T' exp {-^dt, 



(3.24) 



(3.25) 



where z' is defined to be positive on the positive part of the Re (z) axis. 

The far field expressions given below will be expressed in terms of a'^(k cos 6). When 
6-^T, a'^(k cos 6) is singular. It is often useful, in this case, to make the substitution (see eqs 
(2.5) and (2.6)). 

a+{k cos 6) = [l + {a/ik sin 6)]a-{k cos 6). (3.26) 

Power series expansions of (T^(k cos 6) = a'^[k cos 6, (a/k)] and a~(k cos 6) = a~[k cos 6, (a/k)] 
in terms of (a/k), \(a/k)\<C0.1, are given in eqs (2.21) and (2.22). 
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c. Far Fields- Problem P,(0<^<7r) 

In this case, the incident wave originates over the sea. The angle ^o between the direc- 
tion of incidence and the positive a:-axis (see figure 4a) is assumed not to vanish; horizontal 
incidence is therefore excluded. 

Region A: \{2kpy sin [(i9-^o)/2]|»l; 7r>^>(9o+A(9o 

h'^ {exp [ikp cos (^— 6>o)] + oxp [ikp cos id+e^)]] 

2aa^{k cos d) sin 6 exp [i(^p— 7r/4)] 



id'^ {k cos ^o) ((*os ^— cos ^o) {(i-\-ik sin ^o) i^irkpy^ 



(3.27) 



The term in braces is the sum of the incident plane wave and the plane wave whicli is 
reflected from the sea. The last term is obtained by a standard saddle point evaluation. 

Region B: |(2ybp)^sin [(i9-6>o)/2]|</fi«l 

/^^^ (exp [ikp cos (^-^o)]+exp [ikp cos (^+^,)]} ^|^^|^ exp [ikp cos (^-^o)] 

, 2 2(7+(^ cos e) sin e (cos (9+cos ^p) exp [^(^p— 7r/4)] 2c[gi(d,d,eo)—gi(do,e,d)] exp [^(^p— 7r/4)] 
"^ a+ik cos So) ia+ik sin 6) (27rtp)^- "^ ?*o-+(i cos do) b, (6, do) (27rtp)^ 

The first, third, and fourth terms in this expression are continuous across B^=6q. On the 
other hand, the second and fifth are discontinuous across 6=6^ because of the presence factors 
l{e-eo) and [A{efio)lhi{e,do)] respectively (see eqs (3.16), (3.17)). If we write 

1 f(/iM' 

erfc[(Zi(Mo)^p)n = 2-^J^, exp(-f2)rff, 

see eq (3.25)), then it is easy to show that the discontinuity in the second term is just com- 
pensated by that in the fifth term. It should also be mentioned that at d=Bo the expression 
for h^ is independent of (T^{k cos ^o) or G~{k cos ^o). 

Region C: |(2fcp)^ sin [(^-^o)/2]|»i, eo+^eo>e>e,+^^, 

h^^exp [ikp cos (^— ^o)]+^ exp [ikp cos (^+^0)] 

_l 2ao-+ {k cos e) sin 6 exp [i{kp—7r/4)] . . 

io-+ {k cos do) (cos (9— cos do) {a+ik sin 6) (2Trkp)^ ' ^ 

The expression within the braces is the sum of the incident plane wave exp [ikp cos (^—^0)] 
and the plane wave exp [ikp cos (^+^0)] which is reflected off the land; r is the reflection coeffi- 
cient appropriate to the land (cf. eqs (3.24) and (2.2)). 

Region D: |(2^p)* sin [(/3,-^)/2]|«l, 6^6,^0 

h^= {exp [ikp cos (^—^0)]+^ exp [ikp cos (^+^0)]} 

l(e—dj,)2a sin ig&(7+(A: cos g^) exp [ikp cos {^1,-6)] 
ik(f^{k cos ^o)(cos /3ft— cos 6^ cos P^ 

2ca~^{k cos e) sin 6{a—ik sin d) exp [^(fcp— 7r/4)] 
io-+ (k cos ^0) (cos ^— cos ^0) {2Tkp)^ 

. 2ck'[g2{d,d,do)-g2(Pb,dA)] exp [^(fcp-7r/4)] 
"^ ia-^{k cos do)b2(dM (27rkp)' 

+i^W^wM^ """"^ [^•(tp-7r/4)] exp {hkp)ll{dM erfc [(Z,(^,^,)^p)^-]. (3.30) 
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This expression remains valid in region C. Note that the last term contains the factors 
exp {ikp) exp (4fcp) = exp [ikp[\ — 2 sin^ [(jS&— ^)/2]]} =exp [ikpco^ (/3^,— ^)]. As ^ passes through 
^h, ili^yPb) changes sign (see eq (3.21)). From this fact, and eq (3.25), one finds that the jump 
in the third term of eq (3.30) is just compensated by the jump in the last term, the remaining 
terms being continuous across 6=6^- 

d. Far Fields— Problem P,—(i9o=0) 

Again the incident wave originates over the sea, but in the present case the direction of 
incidence is horizontal. 

Region A': |(2A:p)^ sin (i9/2)|»l, 7r>(9>6>,+A^, 

1 79 /•/ Ml aa+(k cos 6) exp [^(/:p+7^/4)] ,_ _^. 
- A^'-'Cxp [ikp cos 6>)+ . ,.j. . .'^.. , .. . .s /o / N^ • (3.31) 

2 ^^ ^ / ' ia+{k) tern {e/2) {a+tk sine) {2wkp)'' 

This expression follows from a straightforward application of the saddle-point method. 
Region B': i(2^p)^sin (^/2)|<i^3<<l and |(2A:p)^ sin [(Pb-0)/2]\<K,<<Cl 

^/.^^exp (*pcos^)+I^^iMO)z:£MO)J exp [^(/:p-7^/4)] 
2 ^ "^""P ^'^^ "^""^^^^ a^(7+(Ar)6l(^,0) (27rA:p)^ 

+^iPwC^ exp [^(Z:p-7^/4)] exp (/,(^,0)/:p)/i(^,0) Xerfc [(l,{dfi)kpy^] 

+ 1(^-^.) ^^^^^J+^^^^^ exp (-i7r/4) exp [ikp cos (/3,-^)] 

F[ff2(^,^,0)-^2(g..^.O)] exM^fe-7r/4)j 
"^ aia+{k)b2(dM {2Tkp)- 

+il^m^ exp [^(A:p-W4)] exp (Z2/:p)/|(^,/3,) Xerfc [(/K^,/?,)^^]. (3.32) 

This expression gives the far field for A* in region A^ also and reduces to the expression given 
in eq (3.3 1) when |(2tp)* sin {d/2)\»l. 

Important simplifications occur when ^=0 — i.e., on the land surface. If the influence of 
the complex pole at ^=Pd is neglected or equivalently if it is assumed that kpo is so large that 

1 (2Z:po)* sin (/3^,/2) |>> 1 then the last three terms are 0[{kpo)~^^^] whenever p>Po. Furthermore, 
the third term just cancels the first term. As a result, it is found that 

Its/ m 2k exp[i(kp—Tr/^)] i\^^^ exp (ikp) . . 

_/, (,^0) — (27rytp)i ^2^2~~7J~- ^^-^^^ 

Here, X is the wavelength in air and c/2 is defined in eq (1.7). This formula was first obtained 
by Griinberg [1942, 1943]. Note that this result is independent of the function a~^. 

On the other hand, when 6=0 and |(2A-p)^ sin (0b/2)\, with p>po, is small compared to 
unity, all the terms of eq (3.32) must be taken into account. The result is 

1,,. ^. 2k r (T+{k COS jg,) cos^ (/3,/2) "] exp [i{kp-ir/4:)] c+jk cos ^,) cos^ (g,/2) 

2 ^ ^^'^ ^ a[_ a+{k) cos 0, J (27ri!:p)* "^ (T+{k) cos /?, 

l+^J^ exp{-ndtj+0[{kp)-"']- (3.34) 

(The terms four and six in eq (3.32) lead, with the help of eq (3.25) to the second term above 
and the remaining terms of eq (3.32) yield the first term above.) Since cos /36=[l+(a/^)^]' 
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and \(afk)\<OA, we have llio following- estimate for h^(p,0) on the land 

1 r 9 r(kp)hl + i) Bin (pf,/2) -| 

2 A^^(p,0) --exp [ip{k'+a')i] |_1 + (^J^ exp (-f^),/^ J+o(a/yt')+0[(/:p)-^/2J. (3.35) 

Tims for sufficiently small \(2kp)^ sin (fib/2)\, h'(p,0) is a gromid wave (cf., eq (1.11)). 
3.3. Some Comments Regarding Coastal Refraction 

As was pointed out by Bazer and Karp [1952], had problem P, been foin minted in terms 
of Ey, the vertical component of the electric field, instead of h (see oqs (1.1)-(1.3)), it would 
have been possible to treat the case of an exciting plane wave incident parallel to the sea but 
making an angle, Jtt— </> say (O<0<7r/2) with the shoreline. (In fact, it was in terms of tlie 
angle </> that G. A. Griinberg [1942, 1943] formulated his discussion of the problem of coastal 
refraction (see section 1).) It would then have been found, emploA'ing tlie method of the 
paper, that Ey varies as (cos <^)'A^(p,0) for \[(k cos 0)p]' sin {I3i,/2)\^y>l on the land, ^ a distance 
p from the shoreline. Here h^ is the expression in eq (3.33). This result is in complete agree- 
ment with Griinberg's and shows, at least within the present theoretical context, that there is 
no refraction at large distances from the shoreline. The possibility of coastal refraction still 
exists, however, when \(kp cos </>)" sin {I3b/2)\<^<il, since, in this case, the presence of a ground 
wave of tlie form exp [^pVF cos^ 0+a^+y!: cos z] would be felt (see eq (3.35)). The attention 
of the authors was drawn to this range by H. G. Booker,^ to whom they are further indebted 
for bringing to their notice an early work by Eckersley [1920] which bears directly on this 
problem. To sum up, the appearance of a ground wave for ''moderate^' values of {k cos </))p 
leads us to believe that coastal refraction may be explained by an extension of the methods 
of our paper. What is required is a vector treatment of the problem for the case of oblique in- 
cidence and a careful investigation of the error terms in the oblique-incidence analogs of the 
results presented in the above summary. 
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